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ABSTRACT 

We use three-dimensional smoothed particle hydrodynamics simulations together with 
a dynamical ray-tracing scheme to investigate the build-up of the first H ll regions 
around massive Population III stars in minihaloes. We trace the highly anisotropic 
breakout of the ionizing radiation into the intergalactic medium, allowing us to predict 
the resulting recombination radiation with greatly increased realism. Our simulations, 
together with Press-Schechter type arguments, allow us to predict the Population III 
contribution to the radio background at ~ 100 MHz via bremsstrahlung and 21-cm 
emission. We find a global bremsstrahlung signal of around 1 mK, and a combined 
21-cm signature which is an order of magnitude larger. Both might be within reach of 
the planned Square Kilometer Array experiment, although detection of the free- free 
emission is only marginal. The imprint of the first stars on the cosmic radio background 
might provide us with one of the few diagnostics to test the otherwise elusive minihalo 
star formation site. 
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1 INTRODUCTION 

One of the most important questions in modern cosmology 
is to understand how the first stars, the so-called Popu- 
lation III ( Pop III), ended the cosmic dark ages at red- 
shifts z < 30 tearkana Loebl l2001 : Bromm k Larson "2004; 
' ICiardi Fer rara 2005^ Their emergence led to a fundamen- 
tal transformation in the early Universe, from its simple ini- 
tial state to one of ever- increasing complexity. The emission 
from the hot, Tgff ^ 10^ K, photospheres of Pop III stars 
began the reionization of primordial hydrogen and helium 
in the intergalactic medium (IGM), although this process 
was co mpleted only lat er on, when more massive galaxies 
formed (|Fan et al.H2006l V In addition, the supernova explo- 
sions that ended the lives of massive Pop III stars distributed 
the first heavy e l ements into the IGM (iBromm et al.ll2003l : 
iGreif et al.|[2007l : iTornatore et~alll2007l : IWise Abelll2008h . 
This latter process might have had a significant impact on 
the physics of early star formation, as metal-enriched gas can 
cool more eff iciently than primordial gas (Bromm & Loeb 
l2003l : lOmukai et al.1 120051 : iJappsen et al.ll2007l . [2009) . 

Based on numerical simulations, a general consensus has 
emerged that the first stars formed in dark matter mini- 
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haloes at z ^ 20 — 30, in isolation or at most as a small 
stellar multiple, and wit h typical masses of ^100 Mq 
(for a recent review, see iBromm et all l2009h . It is crucial 
to observationally test this key prediction. However, it has 
become evident that this will be very challenging. Even 
the exquisite near-IR (~ njy) sensitivity of the upcoming 
James Webb Space Telescope (JWST) will no t suffice to di- 
rectly image such ma ssive, single Pop III stars (|Bromm et al.l 
I2OO1I : [Gardner et"al.l 12006). unles s they explode as ener- 
getic pair- instability su pernovae (|Heger WooslevI I2OO2I : 
IScannapieco et aDl2005l ) . The direct spectroscopic detection 
of recombination line emission from the H 11 region surround- 
ing the Pop III star, as well as from the relic H 11 region 
left behind once the star had died, is beyond the capabil- 
ity of JWST as well, although such line emission might be 
detectable from prim ordial stellar populations inside more 
massive host haloes (jSchaererl I200I l2003l : Ijohnson et al.l 
[2009). 

An alternative approach is to search for the global sig- 
nature from many Pop III stars that formed in minihaloes 
over large cosmic volumes (|Haiman Loeblll997f ) . One such 
probe is the optical depth to Thomson scattering of cos- 
mic microwave background (CMB) photons off free elec- 
trons along the line of sight, determined by the five-year 
Wilkinson Microwave Anisotro py Probe fWMAP) m easure- 
ment to be Te ^ 0.09 =b 0.02 (|Komatsu et al.ll2009l ). This 
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signal, however, is dominated by ionizing sources that must 
have formed closer to the end of reionization, with only 
a smal l contribution from Pop III stars formed i n mini- 
haloes (jGrei f k Bromml l2006l : ISchleicher et al. 2QQ8V A sec- 
ond empirical signature is the combined bremsstrahlung 
emission from the H ii regions in their active and relic 
states around those minihaloes that hosted Pop III stars. 
The resulting free-free radio emission leads to spectral dis- 
tortions that might be detectable in the Ray leigh- Jeans part 
of the CMB spectrum. Recently, the ARCADE 2 experi- 
ment has attempted to measure such a free-free contribu - 
tion from the epoch of the first stars (Ko gut et al.|[2006 lV 
The surprisingly strong signal found, however, cannot orig- 
inate in early Pop III stars, and in any case would over- 
whelm t he much weaker cont ribution from the first stars and 
galaxies (jSeiffert et al.ll2009l ). The most promising detection 
strategy might be to scrutinize the backgroun d from the 
redshifted 21-cm line of neutral hydrogen (Furl anetto et all 
I2OO6I V Once the central Pop III star has died, the relic H 11 
region lef t behind would provi de a bright source of 21-cm 
emission (jTokutani et aLlboOOl l. A gain, individual sources 
are much too weak to leave a detectable imprint, but the 
planned Square Kilometer A rray (SKA) might be able t o 
detect the cumulative signal (|Furlanettoll2Q06l : iLaziolfioOSI V 
We here carry out radiation hydrodynamics simulations 
of the evolution of H 11 regions around massive Pop III 
stars in minihaloes, giving us a detailed understanding of 
the properties of individual sources. We combine this with 
an approximate. Press- Schechter type analysis of the cosmo- 
logical number density of minihaloes as a function of redshift 
to derive the observational signature of the first H 11 regions 
as well as relic H 11 regions, where we specifically focus on the 
free- free and 21-cm probes. We note that we do not include 
the feedback effects exerted by black holes or supernovae, 
which a re possible end products of massive Pop III stars 
(|Heger Woosley 2002). In this sense, we organize our work 
as follows. In Section 2, we describe the simulation setup and 
our implementation of the radiative transfer scheme in the 
smoothed parti cle hydrodynamics (SPH) code GADGET-2 
(Springellliooi). In Section 3, we discuss the properties of 
the first H II regions in their active as well as relic states 
and their observational signature in terms of recombination 
radiation, bremsstrahlung and 21-cm emission. Finally, in 
Section 4 we summarize our results and assess their impli- 
cations. For consistency, all quoted distances are physical, 
unless noted otherwise. 



2 NUMERICAL METHODOLOGY 

Our treatment of ionizing and photodissociating radia- 
tion emitted by massive Pop III star s is very similar to 
the methodol ogy introduced in Johnso n et al.l (2007) and 
lYoshida et aEI (j2007|), with the exception that we here take 
the hydrodynamical response into account, self-consistently 
coupled to the chemical and thermal evolution of the gas. 
This allows us to model dense (D-type) as well as rarefied 
(R-type) ionization fronts, which is crucial for a proper treat- 
ment of the breakout of ionizing radiation. In the following, 
we describe our simulation setup, as well as the numerical 
implementation of the ray-tracing algorithm. 



2.1 Simulation set-up 

We perform our simulations in a cosmological box with lin- 
ear size 200 kpc (comoving), and 256"^ particles per species, 
corresponding to a particle mass of 17 Mq for dark mat- 
ter and c:^ 3 M0 for gas. The simulations are initialized at 
z = 99 with a fluctuation power spectrum determined by a 
A cold dark matter (ACDM) cosmology with matter density 
Qrn = 1 — ^^A = 0.27, baryon density = 0.046, Hubble pa- 
rameter h — ifo/100 km s~^ Mpc~^ — 0.7, where Ho is the 
Hubble expansion rat e today, and spectral index Us = 0.96 
(|Komatsu et al.l [20091 ). We use an artificially high fluctua- 
tion power of (78 = 1.6 to accelerate structure formation in 
our relatively small box, although the cosmological mean 
is given by ag = 0.81. We take the chemical evolution of 
the gas into account by following the abundances of H, H^, 
H~, H2, HJ, He, He^, He^^, and e~, as well as the three 
deuterium species D, D^, and HD. We include all relevant 
cooling mechanisms, i.e. H and He collisional ionization, ex- 
citation and recombination cooling, bremsstrahlung, inverse 
Compto n cooling, and collisional excitation cooling via H2 
and HD (| Glover JappsenllioOTl ). We explicitly include H2 
cooling via collisions with protons and electrons, which is 
important for t he c hemical and thermal evolution of relic 
H II region gas (^over & Abel 2008). 

We run the simulations until the first minihalo in the 
box has collapsed to a density of nn = 10^ cm~^, at which 
point the gas has cooled t o ^ 200 K and beco mes Jeans- 
unstable ( Abe l et al.1 12OO2I : iBromm et alll2002l ). The first 
halo that fulfils this criterion collapses at ^ 20 and has 
a virial mass of ~ 9.4 x 10^ Mq and a virial radius of 
^ 90 pc. Highly resolved simulations have shown that at 
later times, the gas condenses further under the influence 
of self-gravity to nu ^ 10^^ cm~^, where it bec omes op- 
ticall y thick and forms a protostellar seed ( Yoshida et al.l 
I2OO8I ). Due to its residual angular momentum, the central 
clump flattens and likely evolves into an accretion disk. In 
our case, we find a flattened structure already at a density of 
nu = 10^ cm~^ (see Fig. 1). Subsequently, the star grows to 
as ma ssive as ^ 100 Mp wit hin its lifetime of a few million 
years (jBromm Loebll2004l ). However, we note that under 
certain conditions the dis k may fragment to form multiple 
objects of smaller masses (|Clark et al. 2008). Unfortunately, 
the details of the accretion phase and the concomitant ra- 
diative feedback are poorly understood, alth ough some an- 
alytic investigations have been carried out (|Tan McKed 
I2OO4I : iMcKee fanll2008l ). Under these circumstances, it 
seems best to initialize the calculation of the H 11 region at 
the onset of the initial Jeans-instability, when the density 
exceeds nn — 10^ cm~^. 



2.2 Ray-tracing scheme 

The procedure used to calculate the Stromgren sphere 
around the star for a given time-step At i s sim ilar to the 
ray-tracing scheme used in Ijohnson et al.l (|2007f ) . We first 
designate an individual SPH particle as the source of ioniz- 
ing radiation and create a spherical grid with typically 10^ 
rays and 500 logarithmically spaced radial bins around the 
source particle. The minimum radius is set by the smooth- 
ing length of the central particle, while the maximum radius 
is chosen appropriately to encompass the entire H II region. 
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Figure 1. Sequential zoom-in on the first star-forming minihalo at ^ 20. Shown is the density-squared weighted average of the 
hydrogen density along the line of sight, just after the formation of the first Jeans-unstable clump in a 2:^ 9.4 x 10^ Mq minihalo. The 
fiattening of the core due to angular momentum conservation during the collapse is marginally visible, with the consequence that ionizing 
radiation from the central source breaks out anisotropically (see Fig. 3). 



This approach may seem crude compared to e xisting meth- 
ods t hat use adaptive grids (e.g. HEALPIX; iGorski et alJ 
I2OO5I ). but the increased angular and radial resolution to- 
wards the center tend to mirror the existing density profile. 
However, one must proceed with care if the ionization front 
encounters dense clumps far from the source, where the res- 
olution may no longer be sufficient. 

In a single, parallel loop, the Cartesian coordinates of 
all particles are converted to spherical coordinates, such that 
their density and chemical abundances are mapped to the 
bins corresponding to their radius, zenith angle and azimuth, 
denoted by r, ^ and 0, respectively. The volume of each par- 
ticle is approximately given by AV ^ /i^, which transforms 
to Ar = /i, — h/r and Acj) — h/{rsmO). If the volume 
element of a particle intersects with the volume element of 
a bin, the particle contributes to the bin proportional to 
the density of the particle squared. This dependency en- 
sures that overdense regions are not missed if the bin size is 
much larger than the smoothing length, which could occur 
far from the source where the grid resolution is poor. Acci- 
dental flash- ionization of minihaloes is thus avoided. Once 
the above steps are complete, it is straightforward to solve 
the ionization front equation along each ray: 



UnVi 



> dn A^ic 



dt 



47r 



i 



aB / rien+r dr 



(1) 



where n denotes the position of the ionization front, Aion the 
number of ionizing photons emitted per second, the case 
B recombination coefficient, and nn, rie and n+ the num- 
ber densities of neutral particles, electrons and positively 
charged ions, respectively. We assume that the recombina- 
tion coefficient remains constant at its value for 10^ K, which 
is roughly the temperature of the H II and He ill region. 

The numbers of H l/He I and He II ionizing photons are 
given by 



Aion 



dzy 



(2) 



where h from now on denotes Planck's constant, a denotes 
Boltzmann's constant, and z^min is the minimum frequency 
corresponding to the ionization threshold of H i and He 11. 
We assume that massive Pop HI stars emit a blackbody 
spectrum Bu (in erg s""*^ cm~^ Hz""*^ sr""*^) with an effective 
temperature Tes = dex(4.922, 4.975, 4.999) K and lumi- 
nosity L* = dex(5.568, 6. 095, 6.574) Lct^ for a 50, 100 and 
200 M© star, respectively (|Schaererll20"Qi ). This yields 



iVio„,Hi/Hei = [2.80,9.14,26.99] x 10^'' s"^ 
and 

iVion.Heii = [0.72,4.14,15.43] x 10"*® s"' . 



(3) 



(4) 



We do not distinguish between the H ii and He ii region, 
which is a good appro ximation for massive Pop HI stars 
(jOsterbrock Ferlandi r2006) . The lifetimes of the stars are 
given byt* = 3.7, 2.7 and 2.2 Myr, respectively. We neglect 
the effects of stellar evolution, which might lead to a decrease 
of the number of ionizing photons emitted at the end of the 
main sequence (Marigo et al.ll200"ll : ISchaerei]|2002l l , although 
recent investigations have shown that rotating Pop HI stars 
remain on bluer evolutionary tracks and this effe ct might not 
be so strong ( Yoon Langejl2005l : IWooslev Heger . 20061 : 
IVazauez et al.ll2007l ). 

To obtain a discretisation of the ionization front equa- 
tion, we replace the integral on the right-hand side of equa- 
tion (1) by a discrete sum: 



(5) 



where Ar^ is the radial extent of bin i, and the sum extends 
from the origin to the position of the ionization front at 
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the end of the current time- step At. The above equation 
describes the advancement of the ionization front due to 
an excess of ionizing photons compared to recombinations. 
Similarly, the left-hand side of equation (1), which models 
the propagation of the ionization front into neutral gas, is 
discretized by 



nnVi 



>dri 



dt At 



(6) 



where the sum now extends from the position of the ioniza- 
tion front at the previous time-step to its position at the end 
of the current time- step. We perform the above steps sepa- 
rately for the H ii and He ill region, since they require dis- 
tinct heating and ionization rates. For the He ill region, we 
replace the quantities Un and n+ in equation (1) with Un = 
/Hell and n+ = fueiii ^h, where fx is the number density 
of species X relative to nu • We adopt a case B recombination 
rate of — 1-3 x 10~^^ cm^ s~^ for He III recombinations 
to He II (jOsterbrock Ferlandl l2006l l . Applying the same 
prescription to the H II region, we find Un = (/hi + fuei) nu 
and n+ = (/hit + fueii) ^h- Similarly, we adopt a case B re- 
combination rate of aB = 2.6 x 10~^^ cm^ s~^ for hydrogen 
and helium recom binations from their first ion ized states to 
the ground state (jOsterbrock Ferlandl [iooil V We initial- 
ize the calculation of the H ii region at the boundary of the 
He III region, since hydrogen and helium are maintained in 
their first ionization states by rec ombinations of He ill to 
He II (|Osterbrock Ferlandl[2Q06h . We note that the exact 
position of the ionization front is not restricted to integer 
multiples of our pre-defined radial bins, but may instead lie 
anywhere in between. For this purpose we adopt a simple 
linear scaling of the number of ionizations and recombina- 
tions as a function of the relative position of the ionization 
front. The most expensive step in terms of computing time is 
the assignment of the density and the chemical abundances 
to the grid, while the ray-tracing itself requires only a neg- 
ligible amount of time. 



2.3 Photoionization and photoheating 

Once the extent of the H ii and He ill region have been 
determined, the SPH particles within these regions are as- 
signed an additional variable that stores their distance from 
the source. This information is then passed to the chemistry 
solver, which determines the ionization and heating rates, 
given by 



and 



r = r F.a. (l - ^) diy 



(7) 



(8) 



where Fjy and denote the incoming specific fiux and ion- 
ization cross section, respectively. For the case of a black- 
body. 



(9) 



where r is the distance from the source. The resulting rates 
are given by 



_ [0.45,1.32,3.69] x 10"^ _i 
(r/pc) 

[0.42, 1.43,4.29] x 10"^ _i 



^ion,HeI 



{r/pcf 
[0.67,3.72,13.57] x 10" 



(10) 

(11) 
(12) 



[0.40,1.28,3.74] x 10"^" _i _3 

Fhi = nm — - — "2 erg s cm , (13) 

(r/pc) 

[0.41,1.57,4.94] X 10"^" _3 

Fnei = riHei — — "2 erg s cm , (14) 

(r/pc) 

[0.72,4.46,17.13] x 10"^^ _i _3 

FHeii = riHeii — — "2 erg s cm (15) 

(r/pc) 

for a 50, 100 and 200 Mq Pop HI star, respectively. These 
are taken into account every time-step, while the ray-tracing 
is performed only every fifth time-step. Since the hydrody- 
namic time-step is generally limited to one-twentieth of the 
sound- crossing time through the kernel, our treatment of 
the coupled evolution of the ionization front and the hydro- 
dynamic shock is roughly correct. The computational cost 
of runs with and without ray-tracing are typically within a 
factor of a few. 



2.4 Photodissociation and photodetachment 

The final ingredient in our algorithm is the inclusion of 
molecule-dissociating radiation. This effect turns out to be 
of only minor importance in the present study, but will ren- 
der our algorithm capable of addressing a general set of early 
Universe applications. Molecular hydrogen is the most im- 
portant coolant in low-temperature, primordial gas, but is 
easily destroyed by radiation in the Lyman- Werner (LW) 
bands between 11.2 and 13.6 eV. The small residual H2 
fraction in the IGM leads to a very small optical depth 
over cosmological distances, such that even a small back- 
ground can have a significant effect (iHaiman et al.l I2OOOI : 
iGlover Brandl I2OOII : Ijohnson et al.1 12007 ^1. In our imule- 
mentation, we do not take self-shielding into account, which 
becomes important for H2 column densities ^10^^ cm~^ 
(jPraine Bertoldil[T996i l . Such a high column density is dif- 
ficult to achieve in minihaloes, and is more likely to occur 
within the virial radius of the first galaxies (|0h HaimanI 
,.2002 ). However, the onset of turbulence in the first galaxies 
likel y leads to a reduction of self-shieldin g via Doppler shift- 
ing (|Wise Abelll2007l : lGreif et al.ll200"8l ) . For this reason we 
treat the photodissociation of H2 in the optically thin limit, 
such that the dissociation rate in a volume limited by causal- 
ity to a radius r = ct* is given by kB.2 = 1-1 x IO^Flw s~^, 
where Flw is the integral of the specific fiux Fu over the LW 
bands, resulting in 



[1.27,3.38,9.07] x 10" 



(16) 



for a 50, 100 and 200 M© Pop HI star, respectively. We 
equate the photodissociation rate of hydrogen deuteride to 
that of molecular hydrogen. In the present work, we do not 
explicitly include photodetachment of H~ and photodisso- 
ciation of HJ, which might be problematic outside of the 
H II region, where molecules survive collisional destruction. 
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However, in t he context discussed here, this caveat is not 
important f see I Johnson et aDl2007l V 



use the results obtained in this section to determine the re- 
combination signature of the first H ii and He ill regions in 
their active as weh as rehc states. 



3 OBSERVATIONAL SIGNATURE 

In the following, we discuss the direct observational signa- 
ture of the first H ii regions and relic H ii regions in terms 
of recombination radiation, as well as their indirect signa- 
ture in terms of a global radio background produced by 
bremsstrahlung and 21-cm emission. 



3.1 Build-up of H ii and He ill region 

The build-up of the first H ii regions by P op HI stars in 
minihaloes was treated in one dimension by Kitavam a et al.l 
(2004) and Whalen et al. (2004), a nd in three dimensions by 
I Alvarez et al.1 (|2006h . Ubel et al.1 (|2007l ) and lYoshida et al] 
(|2007l ). The latter also treated the build-up of a smaller 
He III region, which is created by the very hard spectrum 
of massive Pop HI stars. The consensus was that recombi- 
nations initially balanced ionizations within the virial ra- 
dius of the host halo, leading to the formation of a D- 
type ionization front. Breakout occured after the density 
dropped sufficiently for the ionization front to race ahead of 
the hydrodynamic shock, becoming R-type. The hydrody- 
namic response of the gas is self-similar, since mini haloes ap- 
proximately resemble sing ular isothermal spheres (|Shu et al.l 
I2OO2I : lAlvarez et al. I I2OO6I ). The relevant parameters are set 
by the temperature of the singular isothermal sphere and 
the H II region, which in our case are T 200 and ^ lO'^ K, 
respectively. 

In F i g. 2, we compare the density profile of the 



et al.l (|2002l ) solution to the simulation for the case of 
a 100 M0 Pop HI star. Interestingly, we find a clear devia- 
tion from the ideal, spherically symmetric solution already 
during the D-type phase, which is caused by the anisotropic 
collapse of the minihalo. Due to angular momentum con- 
servation, the gas spins up and forms a flattened, disk-like 
structure at a density of 10^ cm~^, which can be seen in 
the right panel of Fig. 1, and in the left panel of Fig. 2, 
where it is evident that the density dispersion is almost an 
order of magnitude within the central ~ 10 pc. In response 
to this anisotropy, which is further amplified by the density- 
squared dependence of recombinations, the ionization front 
first breaks out perpendicular to the disk, where the column 
density is lowest. This is visible in the left and middle panels 
of Fig. 3, as weh as in Fig. 2, where the IShu et aD (|2002l ) 
solution is approximately reproduced perpendicular to the 
disk, while the plane of the disk remains neutral and dense. 
Once the ionization front becomes R-type, spherical symme- 
try is asymptotically restored and the H 11 region expands 
to rmi ^ 1.9, 2.7 and 3.7 kpc for the 50, 100 and 200 M© 
Pop HI star, respectively. We find that He 11 ionizing pho- 
tons within the He ill region increase the central tempera- 
ture by a factor of ^ 1.5, leaving only a sm all imprint on 
the d ynamical evolution of the H 11 region f see lYoshida et al] 
[2OO7I). However, the He 11 A 1640 recombination line within 
the He ill region may be used as a distinct probe for the pres- 
ence of massive Pop HI stars (iBromm et"al]l200ll : IOhll200ll : 
iTumlinson et al" ] |200ll : ISchaererll2002l ). In the following, we 



3.2 Recombination radiation from individual H 11 
and He ill regions 

The strongest direct signature of the first H 11 and He ill 
regions is likely generated by recombination radiation, since 
ionizing photons are absorbed by dense gas in the host halo. 
We here concentrate on the Ha and He 11 A 1640 lines, since 
more energetic photons photons are scattered out of reso- 
nance by the neutral IGM, creating extended halo es around 
high-redshift sources (e.g. iLoeb k Rvb"ickil Il999l ). The re- 
sulting fluxes may then be compared to the expected sensi- 
tivity of the Mid-Infr ared Instrument (M IRI) on JWST at 
^ 10 fim wavelengths ([Gardner et al.ll2006l ). The spatial res- 
olution is limited by diffraction, such that a scale of 1 kpc 
at 2; = 20 is marginally resolved, which allows us to approx- 
imate the region of emission as a point source. Using the 
simulation output, the total luminosities are given by 

l^Yloi — jHoi > Je,ijmi,i 

^ pi \ rriH J 

i 



(17) 



and 

^1640 = il640 



I Je,ih 

pi \ mu I 



(18) 



where and ji64o are the emissivity of the Ha and He 11 
A1640 lines at 10"^ K (lOsterbrock h Ferland 2006), X = 0.76 
is the primordial mass fraction of hydrogen, mn is the mass 
of the hydrogen atom, vrii and pi are the mass and density 
of particle i, respectively, and the sum is over all particles in 
the simulation box. From the total luminosity, we determine 
the observed flux with the inverse- square law 



(19) 



where Dl is the cosmological luminosity distance. In Fig. 4, 
we show the observed flux for a 50, 100 and 200 M© Pop HI 
star as a function of time. The emission peaks before break- 
out, when the density in the host halo is still high, reaching a 
maximum flux of 10~^^ erg s~^ cm~^. Once the star turns 
off, the Ha emission drops quite rapidly over the course of a 
few 10 Myr, while the He 11 A1640 emission drops almost in- 
stantaneously in the 100 and 200 cases, due to the high 
recombination coefficient of He ill to He 11. It is noteworthy 
that the emission in the He 11 A 1640 line is generally not 
much lower than that in the Ha line, which m ay be used as 
an indicator for massive Pop H I stars (,Bromm et al. I I2OOII : 
IOhll200ll : ITumlinson et al ] |200ll : ISchaerei]l2002l ). 

For a 10 cr detection with an exposure time of 100 hours, 
the spectrograph on MIRI exhibits a typical limiting sensi- 
tivity of ^ 10"^^ erg s"^ cm"2 (|Panagial I2OO5I ). implying 
that the first H 11 regions are typically five orders of magni- 
tude too faint for a direct detection. We must therefore re- 
sort to indirect methods that rely on their cumulative signal. 
One such signature is the cosmic infrared background (CIB), 
where the redshifted Lya recombination phot ons from z ^ 
10 - 20 might contribute at a detectable level (|Santos et al.l 
I2OO2I : iKashlinskv et al.ll2005l ). Minihaloes, however, are not 
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Figure 2. The hydrodynamic response of the gas to photoheating by a 100 Mq Pop III star after 50 kyr, 300 k yr, and 2.7 Myr (from left 
to right). Shown is the hydrogen density as a function of radius for the simulation (black dots) and the analytic S hu et aP |2002) solution 
(green solid line), as well as the initial density profile of a singular isothermal sphere with 200 K (red dotted line). The functional form 
of the analytic solution is reproduced perpendicular to the disk, where the ionization front breaks out after only a few 10 kyr. However, 
this is not the case in the plane of the disk, where the gas remains neutral and dense until the end of the star's lifetime. 



expected to be important sources for the CIB, as opposed to 
more massive dark matter haloes that host the first galaxies 
([Greif Bromm.2006. ). This leads us to consider the radio 
background as a key diagnostic of the Pop III minihalo for- 
mation site. 



3.3 Radio background produced by 
bremsstrahlung 

The first H ii regions in their active as well as relic states 
also emit bremsstrahlung via thermal motions of electrons 
in an ionized medium. In line with our conlusions of the 
previous section, the signature from an individual source 
is much too faint to be detected. However, the cumulative 
radio signal might be strong enough to be detected by the 
upcoming SKA. We will here fur ther explore this possib ility 
(for a review of earlier work, see iFurlanetto et al.ll2006h . 

Solving the cosmological radiative transfer equation, it 
is straightforward to derive a simple expression for the ob- 
served radio background Ju (in erg s~^ cm~^ Hz~^ sr~^): 

cdt 



Jr 



.0 OT^^"^' ^''^ 

where ti/,o is the present Hubble time and > is the specific 
emissivity of bremsstrahlung, given by 



') (T/10^ 



K 



-1/2 



(21) 



where eff 



10" 



erg s 



cm"^ Hz 



^rie) is the 

volum e- averaged electron density , and T is the temper- 
ature (jRvbicki LightmanI IiQTqI i. We universally assume 
T = 10^ K, since the relic H II region cools quite rapidly 
to - 10^ K via inverse Compton loss es and adiabatic 
expansion once the star has died (e.g. iGreif et al.l l2007l : 
lYoshida et aDl2007l l. Furthermore, we assume i^, = at 
z < 6, sin ce photoheating dur ing reionization evaporates 

to: 



z < b, sm ce pnotoneatmg dur mg reiomzati 
minihaloes (jPiikstra et al.il2004l ). This leads 

r {nl) 

J oo 



ceff 



(1 + 2) 



(22) 



where we relate 
according to: 



3^ to the number density of minihaloes 



^rie) ^ tree nH,b ^HII 



diVps 










dt 



(23) 



Here, tree = (ttBnH,b)~ denotes the recombination time 
for hydrogen atoms, the case B recombination rate for 
T — 10^ K, nH,b the background density, A^ps the num- 
ber of minihaloes per comoving volume, Vhii = A^ion/'n.H,b,o 
the comoving volume of an individual H ii region in its 
active as well as relic state, which is independent of red- 
shift, and A^ion = A^ion t* the total number of ionizing pho- 
tons emitted per Pop HI star (see Section 2). In the above 
equation, we have implicitly assumed that (relic) H ii re- 
gions survive for a recombination time, and that all ionizing 
photons escape into the IGM, which is a good ap proxima- 
tion for massive Pop HI stars in minihaloes (Alvar ez et al.l 
l2006l ). We note that in the range of redshifts considered 
here, the recombination time is larger than the stellar life- 
time and smaller than the age of the Universe. In principle, 
one must also account for the clustering of minihaloes (bias- 
ing) , which reduces the net volume filling factor of H ii re- 
gions IMo k Whitelll996l: llli ev et al. 20031 : iGao et al.ll2005l : 
iReed et al.ll2005l : iGreif Bromm. 20061 ). However, it is ex- 
tremely difficult to determine the importance of this effect, 
since (i) the actual overlap depends on the relative sepa- 
ration of minihaloes, and (ii) previous ionization allows a 
nearby H ii region to become larger than usual. We there- 
fore neglect biasing, but keep in mind that the actual signal 
may be somewhat lower. 

In equation (23), the number density of minihaloes is 
given by 



ATps {z) 



f 



rips (z, M) dM 



(24) 



where rips is the well-kno wn Press- Schechter mass function 
(iPress Schechtei]ll974l ). The minimum mass required for 
efficient cooling wit hin a Hubble time may be found in 
lYoshida et aP (|2003h and lTrenti k Stiavellil (l2009l ): 
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Figure 3. The H ll region created by a 50, 100 and 200 Pop III star (from top to bottom) after 50 kyr, 300 kyr, and at the end of 
their hfetime (from left to right). Shown is the density-squared weighted average of the temperature along the line of sight. The spiral 
structure of the central clump as well as the resulting anisotropic breakout of the ionization front are clearly visible. For increasing stellar 
mass and ionizing photon output, breakout occurs earlier and is more isotropic. Once the ionization front becomes R-type, spherical 
symmetry is asymptotically restored and the H ll region expands to a final radius of rnii — 1-9, 2.7 and 3.7 kpc after 3.7, 2.7 and 2.2 Myr, 
respectively. 
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Figure 4. The observed recombination flux in Ha and He ll A1640 
for a 50, 100 and 200 Mq Pop III star, shown as a function of time 
after reaching the zero-age main-sequence (the black dots denote 
the end of their hfetimes). The emission peaks before breakout, 
when the density in the host halo is still high, reaching a maxi- 
mum flux of ^ 10 erg s"-*^ cm~^. Once the star turns off, the 
Hq; emission drops quite rapidly over the course of a few 10 Myr, 
while the He ll A1640 emission drops almost instantaneously in 
the 100 and 200 Mq cases, due to the high recombination coef- 
flcient of He lll to He ll. The emission in the He ll A1640 line is 
generally not much lower than that in the Ho; line, which is char- 
acteristic for a top-heavy IMF and may be used as an indicator 
for massive Pop HI stars. For a 10 cr detection and an exposure 
time of 100 hours, the limiting sensitivity of the MIRI spectro- 
graph on JWST is approximately 10"-*^^ erg s"-*^ cm~^, indicating 
that the flrst H ll regions are typically flve orders of magnitude 
too faint to be detected directly. 



(25) 



while the maximum mass is set by the requirement that 
cooling must be dominated by molecular hydrogen, i.e. the 
virial temperature must not excee d T 10^ K for atom ic 
hydrogen cooling, resulting in (e.g. iBarkana k Loebll200lh 



Mnaax ^ 2.5 X 10^ Mq ^ ^ ^ 



-3/2 



(26) 



We have found that our results are only marginally affected 
by the upper mass limit, but depend sensitively on the lower 
mass limit, since most minihaloes reside at the lower end of 
the halo distribution function. 

After combining the above equations, we obtain 

C Cff A^ion 



J. 



(27) 



which, for an IMF consisting solely of 100 Mq Pop III stars, 
yields 

J^, 300 mJy sr"^ . (28) 
The brightness temperature, Tb = Ju/2kBi^^ , is given by 



(29) 



In the following, we investigate whether a signal of this mag- 
nitude is observable by the upcoming SKA. 

The sensitivity of radio instruments is generally de- 
fined by the ratio of the effective collecting area Ae to the 
system temperature Tsys- For the SKA with its proposed 
aperture array configuration at low frequencies, Ae/Tsys — 
5 X 10^ m^ K~^ at 100 MHz Q- In this range, the system 
temperature is dominated by Galactic synchrotron emis- 
sion, for which a useful approximation is ^iv en by Tgky — 
180 K (Z//180 MHz)"^-^ (iFurlanetto et al.l l20Q6). resulting in 
Tsvs ^ 800 K and ^ 4 X 10^ m^ The minimum angular 
resolution for an array filling factor of unity at 100 MHz is 
approximately 15 arcmin. At higher resolutions, the sensi- 
tivity decreases much too rapidly for effective imaging. In 
Fig. 5, we compare the sensitivity of the SKA for a 10 a de- 
tection, a bandwidth of Az^obs = 1 MHz, and an integration 
time of 1000 h to the brightness temperature and specific 
flux expected for free-free emission. Although the figure im- 
plies that the free- free signal is detectable by the SKA, we 
have neglected biasing as well as radiative feedback in the 
form of a global LW b ackground, which attenuate s star for- 
mation in minihaloes (| Johnson et al.ll2007l . [20081 ). Another 
complicating issue is the overlap with 21-cm emission, which 
makes it nearly impossible to isolate the contribution from 
bremsstrahlung. In consequence, we do not believe that this 
signal will be observable in the near future. 



3.4 Radio background produced by 21-cm 
emission 

Perhaps the most promising observational signature comes 
from 21-cm emission of the relic H ii region gas once the 
star has died, a prospect that was already investigated by 
Tokutani et al. (2009). An emission signal requires the spin 
temperature Ts of neutral hydrogen to be greater than the 
temperature of the CMB, with its relative brightness de- 
termined by Ts and the size of the relic H ii region. The 
spin temperature is set by collisional coupling with neutral 
hydrogen atoms, protons and electrons, as well as radiative 
coupling to the CMB. Furthermore, it may be modified by 
the so-called Wouthuysen-Field effect, which describes the 
mixing of spin states due to the absorption and re-emission 
of Lya photons (Wouthuysen 1952; Field 1959). The color 
temperature of the Lya background is determined by the 
ratio of excitations to de-excitations, which approaches the 
kinetic gas temperature at high redshifts, where the opti- 
cal depth to Lya scattering is very large (Furlanetto e t al.l 
2006). In this case, adopting the Ray leigh- Jeans approxima- 
tion and assuming Ts ^ T*, where T* = hv2i/kB — 68 mK 
is the temperature associated with th e 21-cm transition, the 
spin temperature may be written as (|Madau et al.|[l997l l 



Ts 



l + yc + ycx 



(30) 



where T^ is the temperature of the CMB. The collisional 
coupUng coefficient yc is approximately given by 



http: / /www. skatelescope.org 



© 0000 RAS, MNRAS 000, 000-000 



The first H ll regions 9 



10 



— — SKA IOct 1000" 
I I I I I I I I I 

200 Mq 

100 Mq 

50 Mq 
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Figure 5. The brightness temperature and specific fiux of the ra- 
dio background produced by bremsstrahlung, shown as a function 
of observed frequency. We have chosen a beam size of 15 arcmin to 
achieve the highest possible resolution and sensitivity at 100 MHz 
for the currently planned configuration of the SKA. The dot- 
dashed, dotted and solid lines correspond to an initial mass func- 
tion consisting solely of 50, 100 and 200 Mq Pop III stars, re- 
spectively. The dashed line shows the sensitivity of the SKA for a 
10 cr detection, a bandwidth of 1 MHz, and an integration time of 
1000 h. Although the free-free signal is in principle detectable by 
the SKA, we have here neglected biasing and radiative feedback, 
which act to reduce the signal. For this reason we do not believe 
that the free-free signal of the first H ll regions in their active or 
relic states will be observable in the near future. 



(31) 



A21T 

where A21 = 2.85 x 10~^^ is the Einstein A-coefRcient for 
the 21-cm transition, and khi and the effective single- 

atom rate coefficients for collisions with neutral hydrogen 
atoms and electrons, respectively. Good functional fits in 
the temperature range 100 K < T < 10^ K are given by 



Km - 10 
and 



-1I/2 3 - 
^ cm s 



Ate = 2 X 10-107^1/2 cm^ s"^ 



(32) 



(33) 



which we have o btained from the rates quoted in 
iKuhlen et~all (|2006[ ). At z<20, the electron fraction in the 
IGM remains above fe =0.1 for most of the lifetime of the 
relic H li region. In this case, the collisional coupling coeffi- 
cient is given by 



103 k; V in / • (34) 

A derivation of the Lya coupUng coefficient i/a requires ra- 
diative transfer of local as well as global Lya radiation, 
which is beyond the scope of this work. We therefore con- 
sider two limiting cases: one in which we only consider colli- 
sional coupling, and the other in which a strong Lya back- 



ground drives the spin temperature towards the gas temper- 
ature (i.e. > 1 or Ts = T). 

The differential brightness temperature with respect to 
the CMB may then be derived as follows. In the Rayleigh- 
Jeans limit and for Ts ^ T*, the monochromatic radiative 
transfer equation for a ray passing through a cloud, evalu- 
ated in its comoving frame, may be written in terms of the 
brightness temperature T^: 



Jo 



where the optical depth at the 21-cm line is given by 



dr 



3c A2inm 
327rzy|^ 



0(z^2i 



Ts 



■ds 



(35) 



(36) 



Here, (/)(z^2i) is the normalized line profile at the resonance 
frequency z/21 and ds is the distance traveled by the ray. In 
our case, the line profile is dominated by thermal broaden- 
ing, with a Doppler width given by 



2A:bT 



(37) 



The amplitude of the line profile at the resonance frequency 
may be replaced by the Doppler width, i.e. (j){h'2i) — Az/^^. 
With this definition, equation (35) yields the differential 
brightness temperature 5Th = — T^, which becomes par- 
ticularly simple for a constant spin temperature and the fact 
that the relic H 11 regions considered here are optically thin: 



5n = (Ts-T^)T 



(38) 



The observed differential brightness temperature is then 
simply given by (5Tb,obs = ^Th/ (1 + z). 

In Fig. 6, we show the observed differential brightness 
temperature for a 100 M© star and the two limiting cases 
discussed above. Note that we have only taken into account 
ionized gas along the line of sight. For collisional coupling, 
the observed differential brightness temperature is of order 
a few 10 mK for ^100 Myr, while for perfect coupling the 
signal is elevated by an order of magnitude to a few 100 mK 
for well over 100 Myr. In reality, the expected signal lies 
between these extremes and is a function of redshift, since 
collisional coupling becomes weaker as the background den- 
sity drops, while Lya coupUng becomes stronger as the Lya 
background rises. At z < 20, where the observationally acces- 
sible signal is produced, the Lya background is likely strong 
enough for the latter to be mo re important ( Fur lanetto 2006, : 
IPritchard FurlanettQll2007l l. 

Next, we determine the radio background produced by 
the integrated 21-cm emission of relic H 11 regions. The dif- 
ferential specific flux observed at the redshifted 21-cm line 
from a single relic H 11 region with differential brightness 
temperature 5Th is given by 

2/cbz^Ii 



(l + z)""" AQ^Tb 



(39) 



where AQ = A/D\ denotes the solid angle subtended 
by the relic H 11 region, A = Trrnn its area, thii = 
(3A^ion/47mH,b)^^^ its radius, and Da the angular diameter 
distance. The average differential specific flux {5Fu) within 
a beam size AQbeam and bandwidth Az^obs is then given by 



dV(z) AzAQbe 



Az/oi 



(40) 
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Figure 6. The observed differential brightness temperature of the relic H ll region around a 100 Mq Pop III star, shown 20, 50 and 
100 Myr after the star has turned off. We delineate the range of possible values by showing the result for collisional coupling only (top 
row), as well as perfect coupling to the Lya background, resulting in z/a ^ 1 or Ts = T (bottom row). In the first case, the observed 
differential brightness temperature is of order a few 10 mK for ~ 100 Myr, while in the second case the signal is of order a few 100 mK 
for well over 100 Myr. We note that the latter is likely more relevant at z < 20, where the observationally accessible signal is produced 
rsee lFurlane'ttQll2006l : IPritchard Furlanettoll200/t ). ^ 



where SF = S Fjy Aiyn / {I + z) , Nps{z) is the Press- 
Schechter mass function defined in equation (24), Az = 
Az^obs (1 + z)'^ /i^2i, and d'^V{z)/dz dQ is the comoving vol- 
ume per unit redshift and solid angle: 



dV(z) _ cDl{l + z) 



dzdQ 



H{z) 



(41) 



where H{z) is the Hubble expansion rate. With the defini- 
tion of the brightness temperature, the average differential 
antenna temperature (^T^) is given by 

{5Th) = — ^ TT( \ r^ii{z) STb(z) . (42) 

Based on our argument above, we assume that the Lya back- 
ground is strong enough for perfect coupling at all redshifts. 
In this case, and for T ^ T^, the average differential antenna 
temperature becomes independent of electron fraction and 
temperature: 



9cM2lT«Afion 



il + z)'/'N^.iz) , 



(43) 



where we have set nm = nu,h in equation (36). We note 
that the observed frequency is related to the redshift via 
J^ohs = 1^21/ {1 -\- z). We have further assumed that the relic 
H II region produced by each star-forming minihalo per- 
sists until the Universe is reionized (i.e. 2; :^ 6), which is 
a good approximation for perfect coupling and T ^ T^. 
Equation (43) thus provides a robust upper limit for the 
collective 21-cm emission from the first relic H 11 regions. 

In Fig. 7, we compare the average differential antenna 
temperature and specific flux for a beam size of A^beam = 
15^ to the sensitivity of the SKA, assuming a 10 a detection, 
a bandwidth of Az/obs = 1 MHz, and an integration time 
of 1000 h. At all frequencies, the maximum 21-cm signal 
from the first relic H 11 regions is of order 10 mK, which 
is well detectable by the SKA. The effects of biasing and 
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Figure 7. The average differential antenna temperature and spe- 
cific flux of the radio background produced by 21-cm emission, 
shown as a function of observed frequency. We have chosen a beam 
size of 15 arcmin to achieve the highest possible resolution and 
sensitivity at 100 MHz for the currently planned configuration of 
the SKA. The dot-dashed, dotted and solid lines correspond to 
an initial mass function consisting solely of 50, 100 and 200 Mq 
Pop III stars, respectively. The dashed line shows the sensitivity 
of the SKA for a 10 cr detection, a bandwidth of 1 MHz, and an 
integration time of 1000 h. In all cases, the 21-cm signal is well 
above the detection threshold of the SKA. The effects of biasing 
and radiative feedback will reduce this signal, but probably not 
enough to fall below the sensitivity of the SKA. Compared to free- 
free emission, the 21-cm signal is typically an order of magnitude 
stronger, and offers the best prospect for indirectly probing the 
first stars. 



radiative feedback will reduce this signal, but probably not 
enough to fall below the sensitivity of the SKA. Compared 
to free- free emission, the 21-cm signal is typically an order 
of magnitude stronger, and offers the best prospect for in- 
directly probing the first stars. Furthermore, the 21-cm sig- 
nal is explicitly frequency-dependent, while this is not the 
case for bremsstrahlung, where a flat spectrum is produced 
(see equation 28). This dependency might allow for a bet- 
ter distinction from other sources of radio emission at these 
wavelengths. 



gions in their active as well as relic states is too faint to be 
directly detectable even with JWST, their collective radio 
emission might be strong enough to be within reach of the 
planned SKA. In particular, we have found that the inte- 
grated free-free emission results in a maximum differential 
antenna temperature of :^ 1 mK, while the 21-cm emission 
is an order of magnitude stronger. Considering the effects 
of biasing and negative radiative feedback, which would act 
to reduce the predicted signal, the free-free signal is likely 
beyond the capability of the SKA, while the 21-cm signal 
will most likely be observable, providing an excellent oppor- 
tunity for indirectly probing the first stars. 

We note that an analysis of the angular fluctuation 
power spectrum will be essenti al to isolate the 21-cm sig- 
nal from other backgrounds (Furlanetto Oh] l2006l l , al- 
though the frequency- dependence of the 21-cm signal might 
already prove useful. Among these are neutral minihaloes, 
which appear in emission due to their enhanced density 
and temperature (|Iliev et al.l |2002|), or IGM gas heated 
by X-rays from super novae (|Ohl I 200 ll ). X -ray binaries 
(iGlover Brandl[20 03) , or the first quasars (|Madau et al.l 
I2OO4I : iKuhlen et aL 12006). A strong absorption signal might 
originate from cold, neutral gas if the Lya background effec- 
tively couples the spin temper ature to the gas temperature 
(jPritchard Furlanetto! l2007l V In addition, there is the sig- 
nal produced by stars (primordial or alread y metal-enriched) 
formed in the first dwarf galaxies (e.g. Naoz & Barkana] 
2008). All of these compete with each other, and more work 
is required to understand their relative importance. One im- 
portant task is to extend the simulations to larger cosmo- 
logical volumes, to measure the aggregate signal from many 
sources in a more robust way. 

Minihaloes may not have been the dominant forma- 
tion sites for primordial stars, in terms of producing the 
bulk of the radiation that drove reionization, or of being 
the source for the majority of the heavy e lements present a t 
high redshifts (jCreif L Broinmll2QQ6l : ISchleicher et al.ll2008l l. 
Nevertheless, they are the ideal laboratory to test our cur- 
rent standard model of the first stars, by providing an ex- 
ceedingly simple environment for the star formation pro- 
cess (Bromm et al. 2009). The next step in the hierachical 
build-up of structure is already highly complex, due to the 
presence of metals, turbulent velocity fields, and possibly 
dynamically signi ficant magnetic fields (IWise Abelll2007l , 
I2OO8I : iGreif et al. I I2OO8I : ISchleicher et al.ll2009l ). It is there- 
fore crucial to empirically probe the minihalo environment, 
and the signature left in the radio background might pro- 
vide us with one of the few avenues to accomplish this in 
the foreseeable future. 



4 SUMMARY AND CONCLUSIONS 

We have introduced a general-purpose radiative transfer 
scheme for cosmological SPH simulations that treats ioniz- 
ing and photodissociating radiation from massive Pop III 
stars in the early Universe. Based on this methodology, 
we have investigated the build-up of the first H 11 regions 
and relic H 11 regions around Pop III stars formed in mini- 
haloes, and predicted their contribution to the extragalactic 
radio background via bremsstrahlung and 21-cm emission. 
Although recombination radiation from individual H 11 re- 
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